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In the modern description of nuclear forces based on chiral effective field theory, four-nucleon 
operators with unknown coupling constants appear. These couplings can be fixed by a fit to the 
low partial waves of neutron-proton scattering. We show that the so determined numerical values 
can be understood on the basis of phenomenological one-boson-exchange models. We also extract 
these values from various modern high accuracy nucleon-nucleon potentials and demonstrate their 
consistency and remarkable agreement with the values in the chiral effective field theory approach. 
This paves the way for estimating the low-energy constants of operators with more nucleon fields 
and/or external probes. 

PACS: 13.20.Gd, 12.39.Fe, 24.80.+y, 24.85.+p, 24.I0.-i, 25.10,+s 



I. INTRODUCTION 



Effective chiral Lagrangians can be used to investigate the dynamics of pion, pion-nucleon as well as nucleon-nucleon 
interactions. In all cases, one has to consider two distinct contributions, namely tree and loop diagrams, which are 
organized according to the underlying power counting JOB]. To a given order, one has to consider all local operators 
constructed from pions, nucleon fields and external sources in harmony with chiral symmetry, Lorentz invariance 
and the pertinent discrete symmetries. Beyond (or even at) leading order in the chiral expansion, these operators 
are accompanied by unknown coupling constants, also called low-energy constants (LECs). In principle, these LECs 
are calculable from QCD but in practice need to be fixed by a fit to some data or using some model [J While m 
certain cases sufficient data exist allowing one to pin down the LECs, often some reliable estimate for these constants 
beyond naive dimensional analysis is needed. In the meson sector, the ten LECs of the chiral Lagrangian at next-to- 
leading order (NLO) have been determined and their values can be understood in terms of masses and coupling 
constants of the lowest meson resonances of vector, axial-vector, scalar and pseudoscalar character, may be with 
the exception of the scalar sector with vacuum quantum numbers |^||. This is called resonance saturation, it has 
been used e.g. to estimate LECs at next-to-next-to-leading order (NNLO) (see e.g. [Q]) or for the extended chiral 
Lagrangian including virtual photons as dynamical degrees of freedom |^] |^. A similar systematic analysis exists for 
the finite dimension two couplings of the pion-nucleon effective Lagrangian [ 11[ , where it was demonstrated the LECs 
are saturated in terms of baryon resonance excitation in the s- and u-channel and t-channel meson resonances. Much 
less is known about dimension three and four couplings, but for certain processes resonance saturation has been shown 
to work quite well, e.g. in neutral pion photoproduction off protons [13]. A somewhat different scheme (including 
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1 These LECs can also be calculated in lattice gauge theory. For a first attempt in the Goldstone boson sector using the strong 
coupling expansion see ref. Q while the most recent quenched calculation for these LECs is given in Q. 

2 For a critical discussion of resonance saturation concerning these LECs, see ||. Note also that the status of resonance 
saturation for the non-leptonic weak LECs is less clear [hoi. 
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also meson-resonance loops) was introduced in the study of the baryon octet masses in ref. |l3). The situation is 
very different concerning few-nucleon systems, where a new type of operators with 2A nucleon fields appears (for 
reactions involving A > 2 nucleons). Only recently, a complete and precise determination of the four S-wave and 
five P-wave (LO and NLO) LECs in neutron-proton scattering has become available flj] , thus it is timely to ask the 
question whether the numerical values of these four-nucleon coupling constants can be understood from some kind of 
resonance saturation.^ This will be the topic of the present paper. As it will turn out, existing one-boson-exchange 
(or more phenomenologically constructed) models of the nucleon-nucleon (NN) force allow to explain the LECs in 
terms of resonance parameters. This paves the way for estimating the low-energy constants of operators with more 
nucleon fields and/or external probes. Before elaborating on the details, we mention that in studies of pion production 
in proton-proton collisions or charge symmetry breaking in the NN interaction, ideas of resonance saturation have 
already been used (T^JI^]. Also, Friar jl^] has discussed aspects of integrating out heavy meson fields to generate 
local four-nucleon operators with given LECs but did not attempt a detailed comparison with existing models of the 
nuclear forces as done here. 

The plan of the article is as follows. In Section || we discuss the effective chiral Lagrangian for nucleon-nucleon 
interactions, in particular the four-nucleon terms and their corresponding couplings constants. We then summarize 
how these LECs are determined at NLO and give a novel prescription to calcuate the NNLO chiral effective field 



theory (EFT) potential. In Section |I| we show how to calculate these LECs from existing boson-exchange or 



phcnomcnological potentials and compare the resulting values with the ones obtained in EFT. Section IV is devoted 
to the study of naturalness of these coupling constants and the implications of Wigner's spin-isospin symmetry. Our 
conclusions are summarized in Section |v|. Some technicalities are relegated to the appendices. 



II. CHIRAL EFFECTIVE FIELD THEORY 



A. Effective Lagrangian and definition of LECs 



To be specific, we briefly discuss the approach to chiral Lagrangians for few-nucleon systems proposed by Weinberg. 
One starts from an effective chiral Lagrangian of pions and nucleons, including in particular local four-nucleon 
interactions which describe the short range part of the nuclear force, symbolically 

C C S — jCtttt + C-kN + £-NN , (1) 

where each of the terms admits an expansion in small momenta and quark (meson) masses. To a given order, one 
has to include all terms consistent with chiral symmetry, parity, charge conjugation and so on. The last term in 
eq.(Q) contains the four-, six-, . . . nucleon terms of interest here. From the effective Lagrangian, one derives the two- 
nucleon potential. This potential is based on (a modified) Weinberg counting more precisely, one organizes the 
unitarily transformed infrared non-singular diagrams according to their power (chiral dimension) in small momenta 
and pion masses (for a detailed discussion, see ref. [fl6||). To leading order (LO), this potential is the sum of one-pion 
exchange (OPE) (with point-like coupling) and of two four-nucleon contact interactions without derivatives. The 
low-energy constants accompanying these terms have to be determined by a fit to some data, like e.g. the two S-wave 
phase shifts in the low-energy region (for np). At next-to-leading order (NLO), one has corrections to the OPE, the 
leading order two-pion exchange graphs and seven dimension two four-nucleon terms with unknown LECs (for the 
np system). Finally, at NNLO, one has further renormalizations of the one- and corrections to the two-pion exchange 
graphs including dimension two pion-nucleon operators. The corresponding LECs can be determined from the chiral 
perturbation theory (CHPT) analysis of pion nucleon scattering. The existence of shallow nuclear bound states (and 
large scattering lengths) forces one to perform an additional nonperturbative resummation. This is done here by 
obtaining the bound and scattering states from the solution of the Lippmann-Schwinger equation. The potential has 
to be understood as regularized, and the regularization is dictated by the EFT approach employed here, i.e. 

V(p,p')^ f R (p)V(p,p')f R { P ') , (2) 



3 Note that in the previous similar work |15| global fits with 26 free parameters where performed, which presumably do not 
allow to pin down the LECs in a unique way. For more details and further discussion on various differences between our 
formalism and the one of the ref. jl5| see Jl^JyJ . 
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where fp(j>) is a regulator function chosen in harmony with the underlying symmetries. Within a certain range of 
cut-off values, the physics should be independent of its precise form and value |pofl . That this is indeed the case has 
been demonstrated in The central object of the study presented here are the LECs related to the four-nucleon 
operators. In a spectroscopic notation these are called Ciso, Ciso, C3S1, C^si, C3D1-351 := Cei, C±pi, C^pq, C3P1, 
and C3P2. In the following, we will collectively denote these as Ci and Q, respectively. The two LECs Ci stem from 
the two momentum independent four-nucleon operators while the seven Ci are related to two-derivative operators as 
they appear in the effective Lagrangian (we have adopted the notation to the two-nucleon potential given in plj ), 

>-NN — L 'NN ' NN + • • • ! 



{ nn = A C s (N^N) (N^N) -I Or (iW, N)^^ N) 



C 

£nn = - \Cx [(N^m) 2 + ((diN^)N) 2 ] - (Ci - ~C 2 ) {N^diN)((diN^)N) 

+ lc 2 (N^N) [N^dfN + dfN^N] 
8 

-|^{[( W 'wr)(WN.)^) + ( W N. W (^ A iV)] 

- {N^N^diN^ajduN) + {N^iN^djN^dkN) 

+ I (i Ca + l Cr ) ^ lk5 ° l + 5a6k ^ + ( 2C3 + \ C ^) SljShl 
x [((didjN^akN) + (N^akdidjN)] (JVV,JV) 

- \ (Ce (fofy + 6 a 6 kj ) + C 4 S l3 5 kl ) (N^kd^^N^mN) 

' f \c 7 {SikSv + 6 a 6 kj ) - (4C 3 - 3C 4 )M«) [(di^a^N) + (d.N^^N)] {N^a t N) , (3) 

where N denotes the (non-relativistic) nucleoli fields, N = (p,n) T , and 07 (/ = 1,2,3) are the Pauli spin matrices, 
and the summation convention for repeated indices is understood. Since we are not considering external sources here, 
we only have partial derivatives acting on the nucleon fields. To arrive at this expression for the most general effective 
Lagrangian with four nucleon field operators, we have made use of partial integration, Fierz transformation and the 
equation of motion for the nucleons. We also require reparametrization invariance p3| of the Lagrangian, which allows 
us to further reduce the number of independent terms as compared to [HJ . The complete derivation of eq. (^) within 
the heavy baryon formalism is presented in |El[ . Note that the effective Lagrangian (|^) corresponds to the rest-frame 
system of the nucleon with the velocity operator given by — (1, 0, 0, 0). The resulting NN contact potential reads 
(in the centre-of-mass (cm.) system): 

Kent = C s + C T {Bl • CJ 2 ) + C X ? + C 2 f + (C 3 q 2 + Cik 2 )^ ■ <7 2 ) + iC 5 ^±^- ■ (qx k) 

+ C 6 {q- ai)(q- a 2 ) + C 7 (k-ai)(k ■ a 2 ) , (4) 

where q = p' — p and k = (p + p')/2 are the transfered and the averaged momentum, respectively, and p (p') 
corresponds to the initial (final) momentum of the nucleons in the cm. system. Closer inspection of eq. ([|) might lead 
to the question why no operators containing the isospin matrices t\ (where i = 1, 2 labels the nucleons) appear? For 
example, /9-meson exchange will naturally lead to a contribution ~ t^-t 2 . In principle, at NLO, one can write down 18 
operators in the effective potential and not just 9 as appear here. What seems to be completely missing are the nine 
operators involving products of isospin matrices. However, we remind the reader that only 9 of these 18 operators are 
independent. The terms in the Lagrangian related to the other 9 can be eliminated using Fierz transformations Jl|,^6[. 
Equivalently, one can perform an antisymmetrization of the two-nucleon potential to eliminate redundant terms as 
used in |23|,[l4|,^lJ . Clearly, the set of operators we choose to work with is one but not the unique possibility. 

As stated before, there are two/seven LECs related to operators with zero/two derivatives. These constants can be 
most easily determined by a fit to the S- and P-wave phase shifts and the 3 S\ — 3 Z?i mixing parameter at low energies, 
which leads naturally to certain linear combinations, i.e. the already enumerated spectroscopic LECs. The precise 
relation of the LECs appearing in the effective Lagrangian to the spectroscopic ones is taken from ref. |l4|] (correcting 
some typographical errors in that reference), 
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Ciso = 47r (Cs — 3Ct) , 

C 150 = 7T (4Ci + C 2 - 12C 3 - 3C 4 - 4C 6 - C 7 ) , 
C351 = 47r (Cs + Ct) , 

C 3 si = I (12Ci + 3C 2 + 12C 3 + 3C 4 + 4C 6 + C 7 ) , 
C 1P1 = y (-4Ci + C 2 + 12C 3 - 3C 4 + 4C 6 - C 7 ) , 
C 3P i = ^ (-4Ci + C 2 - 4C 3 + C 4 + 2C 5 - 8C 6 + 2C 7 ) 
C* 3 p 2 = y (-4Ci + C 2 - 4C 3 + C 4 - 2C 5 ) , 
C 3 po = y (-4d + C 2 - 4C 3 + C 4 + 4C 5 + 12C 6 - 3C 7 ) , 

C 3D1 _ 3S i = C el = (4C 6 + C 7 ) ■ (5) 



B. LECs at next— to— leading order 

Let us now discuss the determination of NLO LECs of the chiral EFT potential. In contrast to what was done 
in 0, we also include the leading charge dependence effect, which is the charged to neutral pion mass difference, 
AM, = M^± — M^o, in the OPE potential (for a systematic study of such effects, see |24j.) Fitting the low neutron- 
proton (np) partial waves (S,P and the triplet S-D mixing) for center-of-mass energies below 50-100 MeV, one 
obtains the numerical values of the LECs for the given regulator and cut-off value. We work here with an exponential 
regulator, 

f R (p) = exp ( V/A 4 ) , (6) 

where the momentum cut-off A is varied between 500 and 600 MeV (a more detailed discussion of various regulator 
functions is given in Jf4|). Therefore, we obtain a range of values for each LEC in the given partial waves. For a 
direct comparison with one-boson-exchange models, we need to further add the TPE contribution which stems from 
the box, triangle and football diagrams. This is done by expanding the contributions of these graphs in terms of 
local operators with increasing powers of derivatives and projecting onto the appropriate partial waves (this method 
is described in more detail below) . The so obtained numerical contributions in each partial wave are listed in table || 
in the column TPE(NLO) (for explicit analytical expressions, see appendix [XJ) . Obviously, these numbers are cut-off 
independent. In contrast, the OPE is retained because all potentials we will compare to include it as well. We note 
that some of these potentials contain a pion-nucleon form factor, but since it only depends on the momentum transfer 
squared and appears quadratically, it does not influence any four-nucleon operator with zero or two derivatives. With 
this in mind, we present in table |the resulting values of the LECs for the cut-off varying from 500 to 600 MeV. This 
is the optimal range found in the study of few-nucleon system fbi]]^ as well as proton-proton scattering in the 
framework used here. Note that in principle we could take smaller values for the cut-off. In such a case one would 
get a slightly less precise description of the data. On the contrary, one could not substantially increase the cut-off 
values if no unphysical deeply bound states are allowed ]f4] , p6] |.p| 

C. Phase equivalent potentials and LECs at next— to— next— to— leading order 

We now turn to the determination of the LECs based on the NNLO potential. Here, we perform a modification as 
compared to the work presented in ref. fl4|| . At that order, the pion-nucleon LECs Ci 3 4 appear, which have been 
taken from the CHPT analysis of ttN scattering in the interior of the Mandelstam triangle |27j. As already shown 
in ref. [fff| , these values can be understood in terms of baryon and meson resonance excitations, with a particularly 
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Although, as will be stressed in section II C, such spurious bound states would not affect low-energy NN observables, the 



direct comparison with the realistic NN potentials would not be possible. 
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strong contribution from the A(1232) resonance. While the natural size for these LECs is 1 GeV -1 , typical values 
found for C3 and C4 from ttN scattering data are C3 = (— 4.70±1.16) GeV -1 and C4 = (3.40±0.04) GeV" 1 , respectively. 
The resulting TPE with insertion of these operators improves the fit but leads to a very strongly attractive central 
potential, as witnessed by the appearance of deeply bound states e.g. in the deuteron channel. These states do, 
however, not influence the low-energy physics in the two-nucleon system. However, the resulting potential is clearly 
not phase-equivalent to the boson-exchange or phenomenological potentials, in which the parameters arc tuned in a 
way that no such additional bound states appear. There is also a more microscopic argument. In a model allowing 
for two-boson exchange (like the Bonn model discussed in refs. j28|) such strongly attractive contributions stem 
from TPE (with intermediate deltas) which are in low angular momentum partial waves almost completely cancelled 
by two-boson-exchange graphs, in which one of the pions is substituted by a p-meson, see e.g. Q and fig. |l|. It 
is even stated in that last reference that "the 27r-contribution appears, in general, too attractive and a consistent 
and quantitative description of all phase shifts can never be reached" . Further work on a detailed understanding of 
correlated np exchange has been performed by Holinde and collaborators, sec p9fl . In the EFT approach, the precise 
order in which the irp diagrams with intermediate deltas start contributing to four-nucleon operators depends on 
the representation of the vector fields (and thus need not appear at the same order as the corresponding 2ir graphs). 
Therefore, we have constructed new chiral potentials at NNLO where in the NNLO TPE graphs we have substituted 
the nN LECs, 

c l ^c i = a - cf , i = 3, 4 , (7) 

using the formalism of ref. p"i"| ] to calculate the c^. More precisely, we have allowed for some fine tuning of the Ci 
within the bounds given in that reference. By this method, the equivalent TPE graphs with intermediate deltas are 
subtracted and the aforementioned cancellations are effectively taken into account. For a typical NNLO fit, we use 
ci = -0.81 GeV" 1 , c 3 = -1.15 GeV" 1 , and c 4 = 1.20 GeV" 1 . From that, we obtain the NNLO TPE contribution 
listed in table || (for explicit analytical expressions, see appendix [X]) . A more detailed description of this procedure and 
further justification of it will be given in a forthcoming publication p0| . The so determined TPE NNLO contribution 
and the corresponding LECs are displayed in table |[ It is important to note that in the cases where the TPE 
contribution is large, the NNLO correction is sizeably smaller than the NLO one. The resulting values for the LECs 
Ci and Ci at NNLO are consistent with the ones found at NLO. That is an important result. 

We are now in the position to confront the LECs determined from chiral effective field theory with the highly 
successful phenomenological/meson models of the nuclear force. Before doing that, some discussion concerning the 
NNLO potential constructed in is in order. It is a perfectly viable scenario to use the unsubtracted values for the 
Ci as done there since the resulting deep bound states do not influence the physics in the two-nucleon system. As noted 
already, a direct comparison of the contact terms in the potential with the ones obtained from the meson-exchange 
or phenomenological approaches can not be made. While the physics in three- or four-nucleon systems does not 
depend on the choice of the unsubtracted or subtracted d , the latter choice is closer to standard nuclear physics in 
which three-body forces lead to small binding energy corrections pp| . In fact, applying directly the potential from 
ref. (IJ] to such systems leads to much smaller binding energies from the two-nucleon forces alone. However, such a 
separation of the total binding energy into NN and 3N contributions is not observable and therefore this scenario is 
not ruled out. These topics will be discussed in much more detail in |3(J. At this stage, both options discussed here 
are viable. It is fair to say that more detailed calculations in few-nucleon systems have to be performed to ultimately 
clarify this issue. We proceed using the modified dimension two pion-nucleon couplings, cf. eq.(R). 



III. LECS FROM BOSON-EXCHANGE AND PHENOMENOLOGICAL NN POTENTIALS 

We consider first genuine one-boson-exchange models of the NN force, in which the long range part of the interaction 
is given by OPE (including in general a pion-nucleon form factor) whereas shorter distance physics is expressed in 
terms of a sum over heavier mesons, 

M=a,p,... 

where some mesons can be linked to real resonances (like e.g. the p-meson) or are parametrizations of certain physical 
effects, e.g. the light scalar-isoscalar cr-meson is needed to supply the intermediate range attraction (but it is not a 
resonance). The corresponding meson-nucleon vertices are given in terms of one (or two) coupling constant (s) and 
corresponding form factor (s), characterized by some cut-off Am- These form factors are needed to regularize the 
potential at small distances (large momenta) but they should not be given a physical interpretation. As depicted in 
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Fig.g, in the limit of large meson masses, keeping the ratio of coupling constant to mass fixed, one can interpret such 
exchange diagrams as a sum of local operators with increasing number of derivatives (momentum insertions). In a 
highly symbolic relativistic notation, this reads, 

(NPiN) {j^T^j iNPjN) = (j^j (NPiN)(NP i N) + (NPN)(NP l N) + ... , (9) 

where the Pi are projectors on the appropriate quantum numbers for a given meson exchange (including also Dirac 
matrices if needed) and Mr is the mass of corresponding heavy meson. It should be kept in mind here that one 
usually makes use of the non-relativistic expansion, i.e. the Dirac spinors on the right-hand side of eq. (^J) coincide 
with the Pauli spinors. In case of a momentum-dependent meson-nucleon coupling, like e.g. for a monopole form 
factor normalized to one at t = 0, 

/(*)=flT3^-T. ( 10 ) 



A 2 - t 



then the coefficient of the first t-dependent term in Eq.(||) is modified to 

9 2 9 2 ( 1 2 



M 4 M 2 \ M 2 A 2 



(11) 



and accordingly for other types of form factors (dipole, monopole normalized to one at t = M R , etc.). The coupling 
constants are either determined in the fit to the NN scattering and bound state data or are taken from other sources, 
the form factor cut-offs always have to be determined from the fit. It is obvious from these considerations that 
such heavy meson exchanges generate four-nucleon terms with zero, two, four, . . . derivatives. In appendix |b], we 
collect the explicit formulae for scalar, pseudoscalar and vector meson exchanges, which can be applied to any of 
the OBE potentials by using the appropriate masses and coupling constants (and should be used instead of the 
symbolic formulae given before) . As a typical example for an OBE potential we consider the Bonn-B variant . Its 
short range part is build from scalar (er, 8), pseudoscalar (rf) and vector meson (ft, uj) exchanges, and the pertinent 
contributions to the LECs Cj, Cj are listed in table ||. Another (more recent) OBE potential is the Nijmegen 93 
potential (denoted Nijm-93) |3^]. The Nijmegen 93 potential is particular since it also includes mesons with strange 
quarks but total strangeness zero (like the scalar e(1300) or the 0(1020)) and a low-energy representation of the 
Pomeron, which usually is needed to describe very high energetic proton-proton scattering. SU(3) flavor sy mm etry is 



imposed so that certain couplings are linked. The various contributions to the LECs are displayed in tabic III. Some 
of the individual terms are unnaturally large (in particular the ones from the Pomeron), but the total contribution of 
the scalar sector is quite similar to the ones in the Bonn-B potential, as comparison of tables f],[II reveals. The vector 



meson contributions (£, uj) are very similar for both potentials. The resulting LECs for these two OBE potentials are 



summarized in table IV. While there is some spread in some of the partial waves, the overall agreement with the LECs 
determined using the chiral EFT potential listed in table [| is rather satisfactory but these results are also somewhat 
surprising, because the phenomenological potential models are not constructed based on any power counting nor chiral 
symmetry, plus in many cases contain quantum field theoretically ill-defined form factors. Still, it is gratifying to 
see that the contact part of the NN potential does not depend on how the short distance physics is parametrized. 
To our knowledge, this is the first time that a direct link between the Weinberg program of systematically deriving 
nuclear forces from chiral Lagrangians to these phenomenologically successful potentials has been achieved in a truely 
quantitative manner. 

There exists also a different class of potentials, which are constructed to give x 2 /datum ~ 1 fits to the NN data 
base like the high-precision charge-dependent CD-Bonn 2000 |33|. It contains two scalar-isoscalar mesons in each 
partial wave up to angular momentum J = 5 with the mass and coupling constant of the second a fine tuned in any 
partial wave. The other high-precision potentials are the Nijmegen I, II |32j as well as the Argonne V18 (AV-18) |3~i| ] 
potentials. For the former, one-pion exchange is supplemented by heavy boson exchanges with adjustable parameters 
which are fitted for all (low) partial waves separately. The AV-18 potential starts from a very general operator 
structure in coordinate space and has fit functions for all these various operators. Note that we have switched off 
the various electromagnetic corrections implemented in the AV-18 potential code. Such type of potentials can also 
be expanded in terms of four-nucleon contact operators with increasing dimension. We do not give the details here 
but only mention that we have done this using numerical methods. The resulting LECs Ci and Cj are also listed 
in table IV. In fig.|| we give a graphical representation of the LECs obtained in EFT, cf. table |, compared with 



the results from the six potential models considered here, see table IV Note that the various LECs are scaled with 



different factors, sometimes including an overall minus sign. This is done to achieve a more uniform representation. 
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Note also that the theoretical uncertainties for the LECs determined in EFT are a) small compared to their average 
values and b) are smaller than the band spanned by the potential models (even if one only includes the high-precision 
ones). 



IV. NATURALNESS OF THE LECS AND WIGNER SYMMETRY 

First, we wish to investigate whether the LECs determined in section || are of natural size. In the present context of 
Weinberg power counting, dimensional scaling arguments allow one to express any term of the effective Lagrangian 
with nucleon and pion fields as well as derivative and pion mass insertions (for a derivation and further discussion, 
see e.g. jl!|) as 

( NH...)N \ l (B*,M*y 2 2 

where the C; m „ are dimensionless numbers and I, to, n are non-negative integers. Here, 21 counts the number of nucleon 
fields, to the number of pions and n the number of derivatives or pion mass insertions. All nucleon isospin operators 
and so on are non-essential to this formula and indicated by the ellipsis. Note that this naive power counting can 
not be applied to cases with spurious bound states, as witnessed by the so-called limit cycle behaviour p5| , |36| , p6[ . 
Here, we only consider potentials with no such spurious bound states, thus the relevant scale for the four-nuclcon 
interactions without derivatives (I = 2, m = n = 0) is the inverse of the pion decay constant, = 92.4 MeV, squared 
and two derivative terms (I = 2, m = 0, n = 2) are suppressed by two inverse powers of the chiral scale A x ~ 1 GeV. 
For the LECs from the Lagrangian (||) naturalness thus amounts to 



Ci~^T, (i = S,T), C^-g™, (* = !,. ..,7), (13) 

Jir Jit ^y 



and the C200 and C202 should be numbers of order one (if there is not some suppression due to some symmetry, see 
below). Such arguments can, of course, not say anything about the signs of the LECs. Also, it is important to realize 
the prefactors which accompany the various terms of the Lagrangian. E.g. there is a relative factor of four in the 
momentum space representation between terms ~ q 2 = (p' — p) 2 and ~ k 2 = (p 1 + p) 2 /A. Such factors need to be 
accounted for. Consequently, we give in table ^ the corresponding coefficients C200 and C202 of the LECs as deduced 
from our NLO and NNLO fits using eqs.(||). Inspection of the table reveals that the numbers fluctuate between 0.3 
and 3.5, i.e. the values found for these LECs are indeed natural, with the notable exception of f 2 Ct, which is much 
smaller than one (except for the upper limit of NLO cut-offs, which is already close to the edge of having stable 
fits, see also the discusion in [^4)). As just mentioned, symmetry can lead to the suppression (or enhancement) of 
certain coupling constants. In fact, 65 years ago Wigner |j7f proposed that SU(4) spin-isospin transformations are 
an approximate symmetry of the strong interactions. Such a transformation has the form 

8N = ie lxu a^T u N , N = ( P \ , /1, v = 0, 1, 2, 3 , (14) 

with = (l,cr), t v — (1,t), and e^ v are infinitesimal group parameters. This symmetry emerges in the large 
number of color limits of QCD |38| and thus features prominently in the nuclear forces derived from Skyrme type 
models. It was recently shown p9| that in the limit where the S-wave scattering lengths a±so & n d a^si go to infinity, 
the leading terms in the EFT for strong NN interactions (with pions treated perturbatively) are invariant under 
Wigner's SU(4) spin-isospin transformations. This can be seen most easily from the leading four-nucleon operators 
as used here, see eq.(||) In this basis, the first term is clearly invariant under Wigner transformations, cf eq.(fl4"l), 
whereas the second term ~ Ct obviously breaks the SU(4) symmetry. In the Weinberg approach employed here, the 
leading order potential consists of these two four-nucleon operators supplemented by the one-pion exchange. Still, the 
Wigner symmetry is kept intact to a good precision since the resulting fit values for Ct are sizeably smaller than the 
corresponding ones for C5 , see table 0. Stated differently, Ct is unnaturally small because of the Wigner symmetry. 
This can be understood from the fact that at very low energies, where one is essentially sensitive to the (S-wave) 
scattering lengths, the pion-exchange contribution can be expanded in powers of momenta, leading to terms with 
at least two derivatives, see appendix |b|. One thus effectively recovers the situation eluded to in ref. |59|. However, 
for larger momenta (say of the order of the pion mass) , the non-perturbative treatment of the pions as proposed by 
Weinberg is mandatory. 
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V. CONCLUSIONS AND SUMMARY 



In this manuscript, we have investigated the low-energy constants with zero and two derivatives that appear in the 
four-nucleon contact interactions of the chiral effective Lagrangian for the nucleon-nucleon forces. Our main findings 
can be summarized as follows: 

(1) We have determined the LECs for the NLO and NNLO potentials, including the dominant charge-dependence 
effect from the pion mass difference in the one-pion exchange. To avoid the unphysical bound states at NNLO, 
we have argued that one has to subtract the delta contribution from the dimension two pion-nucleon LECs. 
This is in agreement with two-boson-exchange models, where the two-pion-exchange contribution is cancelled 
largely by up graphs. 

(2) We have shown how to deduce similar type of contact operators from boson-exchange models in the limit of 
large meson masses. This allows to calculate the LECs in terms of meson-nucleon coupling constants, meson 
masses and (unobservable) cut-off masses. In a similar manner, one can examine the so-called high-precision 
potential models. We have found that in all cases, the LECs determined from these models are close to the 
values found in EFT, which can be considered as a new form of resonance saturation. 

(3) We have shown that with the exception of one dimension zero coupling (the LEC Ct), all LECs are of natural 
size. The smallness of Ct is due to Wigner's spin-isospin symmetry, as was already pointed out for the case of 
a theory with pions integrated out or treated perturbatively. 

Clearly, these findings have further-reaching consequences. On one side, they might allow to further constrain models 
of the nucleoli nucleoli interaction applicable at energies where the EFT description can not be used. On the other 
hand, in case of external sources (like e.g. photons) or multi-nucleon operators (as they appear e.g. in the description 
of the three-body forces), these considerations will allow to at least estimate novel LECs that will appear. In the latter 
case of three- and more-nucleon systems performing a direct fit with 5 adjustable parameters (if the leading non- 
vanishing three-nucleon force is included) to 3N observables will be a very expensive task with respect to computer 
power. Therefore it might be very helpful to have a rough estimation for the values of various couplings appearing in 
the 3N force. 
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APPENDIX A: REDUCTION OF THE TWO-PION EXCHANGE CONTRIBUTIONS 

As stated before, we have to add the contribution of the TPE to the LECs so as to be able to compare with the boson 
exchange potentials. The explicit expressions for the renormalized TPE potential at NLO can be found in ref. p6| . 
Expanding those in powers of q and k allows for a mapping on the spectroscopic LECs (of course, the TPE contains 
many other contributions, which are, however, of no relevance for this discussion). We get 

_1 ~ NLO = {l + Ag\~Sg\)Ml 

3 3S1 24tt/4 
2 + \lg\ - 88g\ 

1447T/4 

2 + 17 g\ - 40g\ 

487T/4 

A 

2 + \7g\ - 16g\ 

727T/4 



ANLO 
^lSO 



/.NLO 



✓-.NLO 
°3S1 



/.NLO 



pNLO _ 
°1P1 



^.NLO 



c 



n 

L ■ 



2 + 17g| + 74g 4 
216vr/4 

XL() _ 2 + 17g| - 61g 4 

2167T/4 



3P1 



NLO 
3P2 



2+17.g2 _7. 9 4 



216*/* (A1) 

Note that in the chiral limit, the two leading contact interactions do not get renormalized by TPE. Furthermore, these 
expressions only depend on the lowest order pion-nucleon coupling ~ g A (or, by virtue of the Goldberger-Treiman 
relation, on g^NN)- Similarly, we can give the additional TPE NNLO contributions to the various LECs (for an 
explicit expression of the renormalized NNLO TPE potential, see e.g. ref. [H), 



/^NNLO 
°1S0 



ANNLO 
°3S1 



A.NNLO 
°1S0 



°3S1 — 



gj(-16 + 192m(-2ci + c 3 ) + 25g A )M% 
256m/ 4 
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/-rNI> 



°1P1 



^NNLO 
°3P0 



r NNLO 
°3P1 



256m/ 4 




5^(-368 - 192m(10ci - llc 3 + 4c 4 ) 
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3072m/ 4 




<?i(16(-7 + 4m(10ci - lies + 4c 4 )) + 81 5 i)M w 


1024m/ 4 




g A (8 + 32mc 4 - 7g A )M* 




64 v / 2m/ 4 




g A (-368 + 64m(10ci - llc 3 - 12c 4 ) 


+ 305 5 i)M w 


1536m/ 4 




g A (176 + 192to(10ci - llc 3 - 8c 4 ) + 


mig 2 A )M„ 


4608m/ 4 




^(464 + 192m(10ci - llc 3 + 2c 4 ) - 


545g 2 A )M„ 



L ■ 



NNLO 
3P2 



4608m/ 4 

ffi (112 + 192m(-10ci + llc 3 + 2c 4 ) + 28\g 2 A )M 7 , 



4608m/ 4 



(A2) 



with m the nucleon mass. These expressions depend on the dimension two LECs ci l3)4 as discussed before. We note 
that all these contributions vanish in the chiral limit. 



APPENDIX B: REDUCTION OF ONE-BOSON-EXCHANGES 



Here, we give the explicit expression for scalar, pseudoscalar and vector meson exchange contributions to four- 
nucleon operators with zero or two derivatives, as depicted in fig. ||. Note that we will also include 1/m as well as 
1/m 2 corrections, which are, strictly speaking, of higher orders in the power counting scheme we are working with 
and not (or partly) present in the NLO and NNLO potentials. There is, however, no contradiction since adding or 
subtracting those terms from the potential would lead to changes smaller than the level of accuracy of our approach. 
The contributions for a particular OBE potential can be obtained by using the appropriate masses, coupling constants 
(and form factors) employed there. To obtain the most general expressions, we include any form factor as 

F M (q 2 ) =a 1+ a 2 ^- + 0(q 4 ) , (Bl) 

where the coefficient a.\ = 1 if the form factor is normalized to one at q 2 — or a 4 ^ 1 if the form factor is 
normalized to one at q 2 = —M^ (with Mm the mass of the meson under consideration) or it might include the 
meson-nucleon coupling constant gu- We will give generic expression where the corresponding vertices are written 
as <7m x Fniq 2 ), with F M expanded as just discussed. Note for the exchange of isovector bosons like the ir or the 
p, the given expressions have to be multiplied by a factor t\ ■ r 2 leading to a factor of (-3) for the T = potential 
considered here. 
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1. Scalar meson exchange 



The Lagrangian for coupling of an scalar-isoscalar meson with mass Ms and coupling constant gs reads 

Cs = gsi>i><t> , 



(B2) 



where %p denotes the relativistic nucleon field and <p the scalar meson (for an isovector, one simply replaces <p by 
t • <fi, with r k (k = 1, 2, 3) the usual Pauli isospin matrices). In the non-relativistic expansion, the momentum space 
expression for the corresponding exchange potential with a form factor (if applicable) characterized by the cut-off As 
reads up to terms of order 1 /m 2 



Vs(q,k)=- 



9% 



Ml 



1 - 



q 



2m 2 8m 2 2m 2 



S-(qxk) 



(B3) 



where S = ((?i + <? 2 )/2 is the total spin of the two-nucleon system. The fully relativistic form of this exchange can 
be found e.g. in ref. |pl|. This gives the following contributions to the spectroscopic LECs: 
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U 3P1 
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3P2 
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4irg s al 
Ml ' 




°3S1 — 
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MJA 2 










2-Kg 2 s ai 


((-4 + M|/m 2 )ai 


f 8M|a 2 /A| ) 




3M| 




2irg s ai 


(8m 2 Mla 2 ~ (4m 2 


- 3M 2 ) ai A 2 ) 




3m 2 Mj A| 




47T5|ai 


((-2 + M 2 /m 2 ) ai 


f 4M|a 2 /A 2 ) 



3MJ 



8-Kgla! (2M|a 2 - aiA|) 



3M| A 2 



(B4) 



2. Pseudoscalar meson exchange 



The Lagrangian for coupling of a scalar-pseudoscalar meson with mass Mp and coupling constant gp reads 

£p = -gptprfipTr , (B5) 

where tt denotes the pseudoscalar meson (for an isovector, one simply replaces it by r ■ tt). This is the so-called 
pseudoscalar coupling. Equivalently, one can also use a derivative type (pseudovector) coupling 



C!p = -^7V^ . (B6) 



At tree level, these couplings are equivalent provided gp/m = fp/Mp. Of course, chiral symmetry enforces the 
derivative coupling for the Goldstone bosons. In the non-relativistic expansion, the momentum space expression for 
the corresponding exchange potential with a form factor (if applicable) characterized by the cut-off A p reads up to 
terms of order 1 / m 2 

^---i-i^W 1 ^ 1 - (B7) 



Again, the fully relativistic form of this exchange can be found e.g. in ref. |31|. This gives the following contributions 
to the spectroscopic LECs: 
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S^YP 
2 1P1 



3P2 



o , 



3P0 



3 n p 

i°3Pl 



(B8) 



with 



r P = 



2 2 

m 2 M 2 p 



(B9) 



3. Vector meson exchange 



The Lagrangian for coupling of a vector meson with mass My and coupling constants gy (vector coupling) and fy 
(tensor coupling) reads 



l v = -gvftfM* - p-^^a^ ~ a^) , 

Am 



(BIO) 



where ^ denotes the isoscalar-vector meson (for an isovector, one simply replaces (j)^ by r-0„). In the non-relativistic 
expansion, the momentum space expression for the corresponding exchange potential with a form factor (if applicable) 
characterized by the cut-off Ay reads up to terms of order 1 / m 2 



V 



2m 

•2 
V 



3fc 2 q 2 3i * _ r , _ _ <f 2 1 _ _ 
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mm mm 



A 2 
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(Bll) 



Again, the fully relativistic form of this exchange can be found e.g. in ref. |31|. This gives the following contributions 
to the spectroscopic LECs: 
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TABLES 
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TABLE I. Values of the LECs at NLO and NNLO for the cut-off values A = 500 . . . 600 MeV. Also given are the contributions 
from two-pion exchange at NLO and NNLO which are contained in the values of the LECs as explained in the text. The C; 
are in 10 4 GeV" 2 and the C t in 10 4 GeV~ 4 . 
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TABLE II. Contributions of the various boson exchanges to the LECs for the Bonn-B potential and the corresponding sum. 
The Ci are in 10 4 GeV~ 2 and the C % in 10 4 GeV~ 4 . 



13 



LEC 


?7 


W 
V 


P 






Q>0 




Jo 


Q>2 




Ciso 


u.uuu 


000 
u.uuu 


n 090 

U.UZU 


U.ZO f 


nni 

U.UUl 


O 0*31 
— U.UOl 


n ^78 


n 901 

— U.ZUl 


nni 

U.UUl 


O /1QO 

u.'iyu 


Ciso 


n o/i 1 

U.U41 


n ni *3 

U.U10 


n 1 qi 
u. iy 1 


n /i /i ^ 


009 
— u.uuz 


U. 10^ 


O.'lOl 


n 8^7 

U.oO i 


nn^ 

— U.UUO 


9 89Q 
— z.ozy 


^3S1 


n nnn 
u.uuu 


n nnn 
u.uuu 


0^ 


n 9^7 


nni 

U.UUl 


O OQ/1 


n ^78 

— U.O/o 


n 901 

— U.ZUl 


00^ 

— U.UUo 


O /1QO 

u.'iyu 


L'SSl 


n m a 


n nn/i 

— U.UU^ 


n ^0 
u.oou 


n 700 

— U. 1 uu 


on*} 

— U.UUO 


/in^ 

— U.^tUO 


0.401 


n 8^7 

U.oD i 


O 01 ^ 
U.U10 


9 89Q 
— z.ozy 


C e i 


-0.038 


-0.012 


-0.383 


0.090 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


Cipi 


-0.027 


-0.009 


0.431 


0.423 


0.002 


0.250 


-2.194 


-0.539 


-0.010 


1.791 


C3P0 


-0.082 


-0.026 


0.645 


1.167 


0.006 


-0.072 


-1.985 


-0.466 


0.003 


1.613 


C3P1 


0.054 


0.017 


0.032 


0.660 


0.003 


-0.078 


-2.087 


-0.501 


0.003 


1.700 


C3P2 


0.000 


0.000 


0.114 


0.457 


0.002 


-0.090 


-2.308 


-0.579 


0.003 


1.887 



TABLE III. Contributions of the various boson exchanges to the LECs for the Nijmegen 93 potential. Pseudoscalars: rj,rj' , 
vectors: p,uj,<f>, scalars: a , e, / , a 2 , Pomeron. The C; are in 10 4 GeV -2 and the Ci in 10 4 GeV -4 . 
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TABLE IV. Results for the LECs for various OBE and other type of potentials as explained in the text. The so-called high 
precision potentials are marked by a star. The d are in 10 4 GeV~ 2 and the d in 10 4 GeV" 4 . 
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TABLE V. Naturalness coefficients of the LECs at NLO and NNLO for the cut-off values A = 500 . . . 600 MeV. The / 2 Cs,t 
and the / 2 A 2 d are dimensionless. 
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FIG. 1. Classes of diagrams which cancel to large extent. One representative TPE graph and one np graph are shown. Solid, 
double, dashed and wiggly lines represent nucleons, deltas, pions and p-mesons, respectively. 




FIG. 2. Expansion of a meson exchange diagram in terms of local four-nucleon operators. The dashed and solid lines denote 
the meson M — p, a, u, ... and the nucleons, respectively. The blob and the square denote insertions with zero and two 
derivatives, in order. The ellipses stands for operators with more derivatives. 
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FIG. 3. Graphical representation to compare the LECs determined in EFT with the values obtained from the various models. 
The left bar refers to NLO, the middle bar to NNLO and the filled circles to the various potentials discussed in the text. Note 
that the units are arbitrary, i.e. different scale factors (including in some cases an overall minus sign) have been assigned to 
the various LECs to obtain a more uniform representation. 
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